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Abstract 

We discuss D-dimensional scalar field interacting with a scale invariant 
random metric which is either a Gaussian field or a square of a Gaussian 
field. The metric depends on d dimensional coordinates (where d < D). 
By a projection to a lower dimensional subspace we obtain a scale invariant 
non- Gaussian model of Euclidean quantum field theory in D — d or d 
dimensions. 
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1 Introduction 



We consider a new method of a construction of Euclidean fields. A scalar field 
in D dimensions is interacting with a metric depending on d dimensional coor- 
dinates. An averaging over the metric and a projection of the scalar field to an 
s dimensional subspace leads to a scalar field which is Euclidean invariant in R s 
(we consider s = D — d and s — d). If the metric field is scale invariant with a 
scaling dimension 2j then the scalar field is also scale invariant with a scaling 
dimension depending on 7. We discuss two models for the random metric. In 
the first model we consider a square of a Gaussian random field. We are unable 
to derive an upper bound for correlation functions in this model. Then, we 
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consider a metric which is Gaussian. We obtain scale invariant correlation func- 
tions with explicit upper and lower bounds. Our primary interest in this class of 
models [0-0 comes from quantum gravity. However, the method may be useful 
for a construction of relativistic quantum fields (although at the moment we are 
unable to prove the crucial Osterwalder-Schrader positivity ||). The model can 
be interesting for statistical physics as a continuum version of spin glass models 
[|J. The lattice version of our model describes spins with a random coupling 
between them which is either Gaussian (then we have a mixing between ferro- 
magnetic and antiferromagnetic couplings) or a square of the Gaussian field. A 
calculation of the average over the random coupling can be done explicitly. As 
a result we obtain models with many interacting spins in contradistinction to 
the conventional models based on bilinear spin-spin interactions. 

2 /^-dimensional scalar fields 

We consider a complex scalar matter field $ in D dimensions interacting with 
gravitons varying only on a <i-dimensional submanifold. We split the coordinates 
as x — (X, x) with x e R d . Without a self- interaction the $<!>* correlation 
function is equal to an average (W(g) is the gravitational action) 



over the gravitational field g of the Green's function A 1 (x,y) of the operator 




(1) 





k=D-d 
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In order to calculate the average (1) we repeat some steps of refs.§-@. We 
represent the Green's function by means of the proper time method 

1*00 

A-\x,y)= dT(exp(TA))(x,y) (3) 



Jo 

For a calculation of (exp (tA)) (x, y) we apply the functional integral 

K T (x,y) = (exp (tA)) (x, y) = J Vx exp(- \ J f f - | / gi»> (x) **f **? ) 
S (x (0) - a;) 5 (x (r) - y) 

(4) 

In the functional integral (4) we make a change of variables (x — > b) determined 
by Stratonovitch stochastic differential equations || 

dx n (s) =e%(x(s))db A (s) (5) 

where for fi = 0, 1, D-d- I 



ande5 = ^ if 0>£>-d-l. 

After such a change of variables the functional integral in eq.(4) becomes 
Gaussian. In fact, this is the standard Wiener integral and b^if) for each A are 
independent Brownian motions 

E[b A (t)b c (s)} = S AC mm(s,t) 

The solution q T of eq.(5) consists of two vectors (Q,q) where 

q(r,x) =x + b(r) (6) 
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and Q has the components (for /x = 0, D — d — 1) 

Q"(t, X) = + f < (q (a, x)) dS a ( S ) (7) 
Jo 

The kernel is 

K T (x,y)=E[6(y-q T (x))] 

= E[5(y - x - b(r)) «J - Q M (r, X))] 

Using eq.(7) and the Fourier representation of the (5-function we write the kernel 

K T in the form 

K T (x,y) = (2n)- D + d fdPexp(iP(Y-Xj) 

E[S (y - x - b (r)) exp (-* / P„e£ (q {a, x)) dB Q (s))] W 

We may choose a Gaussian field as a model for the tetrad (as we did in ref. [jy ) 

(e£(xK(y)> = T^(x - y) = o£|x - y^ (9) 

where a is a scale invariant tensor. Then 
(iT r (s, y)) = (2ir)- D+d J dP exp (»P (Y - X)) 

E[S (y - x - Vrb (1)) exp (-r^P, dB« (s) f ° dB< (s') T% (b (s) - b (s')) 

(10) 

where we have changed the time s — » ts , used the equivalence b(rs) = y/rb(s) 
and the scale invariant form of the two-point function (9). Moreover, we renor- 
malized the kernel K T removing from it the term (see 00) 



exp(--rr^(0)P / ,P,) 



It can be seen that this procedure is equivalent to the normal ordering of the 
metric as a square of the tetrad 

<T{x) = <(*)<(*) ->: e£(x)e v a (x) := e£{x)e? a (x) - (e£(x)e» a (x)) (11) 



We can prove that the double stochastic integral in eq.(10) is a finite square 
integrable random variable if 27 < 1. However, it remains unclear whether the 
momentum integral in eq.(10) is finite. 

We can work without the stochastic integrals (8) if we explicitly integrate 
over B. The random variables b and B a are independent. Hence, using the 
formula || 

£[expi j f a (q)dB a ] = £[exp(-i / f a f a ds)] 
we can rewrite eq.(8) solely in terms of the metric tensor 



K T (x, y) = (2Tr)- D+d J dP exp (iP (Y - X)) 

E[5 (y - x - b (t)) exp (-± f (q (s, x)) P v ds)] 

Let be the characteristic function of g^ v 



(12) 



J(h) = {exp(-~g(h)\) (13) 

Then, the mean value of the kernel (12) can be expressed in the form 

(K T (x, y)) = (2n)- D + d J dP exp (iP (Y - X)) 
E[5(y-x-b(T))J(h)} [M) 

where g(h) = J dzg^" (z)/i Ml/ (z) and 

h^(z)=P^ l P l/ [ 5(z-x-b(s))ds 
Jo 

If is Gaussian then J can be calculated explicitly 

(e X pf-i<7(/A=det(l + /ir)-* (15) 

where on the r.h.s. the renormalization of the determinant (through a multi- 
plication by exp (iTr (r/i)) defining det2, see 0) is equivalent to the normal 
ordering (11) (and subsequently to the renormalization of the kernel (10)). 



We consider next another (simpler) model where the metric is Gaussian with 
two-point correlations 

(<r W(y)) = -D»»w( x - y) = -C™|y - ^ (16) 

where C (a scale invariant operator) must be positive definite if the momen- 
tum integrals in the final formula are to exist. This requirement is not satis- 
fied in a linearized Einstein gravity || ( e.g., in the transverse-traceless gauge 
PnPr g ni (x) would be zero in a covariant Z?-dimensional gravity; however our 
gravity is d-dimensional) . The conformally flat metric C' AW;<Tp = 5^ l, 8 ap would 
be a satisfactory model for our purposes . 

The average over g in eq.(12) can be calculated 

(K T (x, y)) = {2n)- D+d J dP exp (»P (Y - X)) 

B[*(y-x-Vrb(l)) (17) 
exp (-ir 2 - 2 ^ f X P^PvPpD^P (b (s) - b («')) dads')] 

(as in eq.(10) we have changed the time s — > rs). By a scaling of momenta we 

can bring the propagator of eq.(3) to the form 

/>oo 

(A- 1 (x,y))= dTT-i- ( - D - d ^ 1 -^/ 2 F 2 {T-i(y- X ),T-i+i{Y - X)) (18) 
Jo 

3 A projection to D — d dimensions 

The two-point function (18) has a different scaling behaviour in x and X di- 
rections. We obtain a fixed scaling behaviour setting x = y = 0. Then, we 
have 

{A- X {x, y)) = R\X - Y |- D+2 -T^(d-2) (19) 
G 



where R is a positive constant. Hence, if all the correlation functions are scale 
invariant then 

*(0,X) ~ A^ + t^t^$(0,AX) (20) 

where the equivalence means that both sides have the same correlation functions. 

In order to prove that R is finite and not zero we need upper and lower 
bounds for the Gaussian model (17). We show first that the bilinear form 
(fj, (A^ 1 )^) is finite and non-zero on a dense set of functions /. For this 
purpose we choose 

/ k (X) = (27ra)-* exp(-^X 2 + zkX) 

Then, (we keep x ^ y in order to show that the model of sec. 2 is non-trivial; 

for a scale invariant model of this section x = y = 0) 

(/k.M-VkO = (2n)- D +«f °°dTT-ifdP 
E[5 fr-5(y-x)-b(l)) 

exp(-i(P-k) 2 -^(P-k') 2 (21) 
_ i T 2-2 7 £ p^p a p uPp D^P (h (s) - b (*') )dsds')] 

In our estimates we apply Jensen inequalities in the form (for real functions A 

and /) 

E[expA]>expE[A] (22) 

and 

£[exp(- J dsds'f{s,s')^]<J^ dsds'E[exp(-f(s,s'))] (23) 
An upper bound can be obtained by means of the Jensen inequality (23) ex- 
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pressed in the form 
(fit, fa) 

< 2 J Q dr f Q ds J Q ds' J duidu 2 dP 

r-i exp (-£ (P - k) 2 - £ (P - k') 2 ) p{s', n 1 )p(s - a', u 2 - u x ) ( 24 ) 
P (l - 8, r-5 (y - x) - u 2 ) exp (-^P^P^i?™ ( Ul _ U2 )) 

where p(s, u) = (27rs)~^ exp(— u 2 /2s). We can convince ourselves by means of 

explicit calculations (using a proper change of variables) that the integral on 

the r.h.s. of eq.(24) is finite. For the lower bound it will be useful to introduce 

the Brownian bridge || starting from and ending in u defined on the time 

interval [0, 1] 

a(u, s) = us + c(s) 

where c is the Gaussian process starting from and ending in with mean 
equal zero and the covariancc 

E[cj(s')c k (s)} = S jk s'(l - s) 

for s' < s. Then, the S function in eq.(21) determines the Brownian bridge and 
the Jensen inequality (22) takes the form 

(/k, {A-i)h,) > (2n)- D + d /~ drr-i f dP 

"3 r2 ~ 27 Jo P^P*PvP p E[D^°p (a(r-3y-r-5 X ,s) - a (r~*y - r^x, A )dsds'} 



(25) 

where the expectation value in the exponential on the r.h.s. of eq.(25) is equal 
to 

JdufdsJ°ds'(2™(s,s'))-* exp^^^u 2 ) (2g) 
|u — r _ 2 S (y — x) + r _ 2 s ' (y — x) |~ 47 
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where lu(s, s') = (s — s')(l - s + s'). It is finite if 7 < \ (the form (16) of the 

graviton two-point function is assumed). 

We compute now higher order correlation functions in the Gaussian model 

<*(a;)*(a; / )**(w)**(w , )> 

= (A- 1 (x,y)A- 1 (x',y')) + {x ^ x') 

= (2n)- 2D+2d J d Tl d T2 J dPdP' exp (iP (Y - X) + iP' (V - X')) 

E[5 (y - x - b (71)) S (y' - x' - b' (r 2 )) 

exp ( - i J T1 J T1 P^PvPpD^P (b (a) - b (a')) dads' 

-3 Jo' 2 Jo' 2 P'^PIP' P D^" (b' (a) - b' (*')) dsds' 

-3 Jo' 1 Jo' 2 P^uP'^D^P (x - x' + b (s) - b' (s')) dads')] + (x - x') 

(27) 

where (x — ► x') means the same expression in which x is exchanged with x' . 
The fourlinear form (27) calculated on the basis / reads 

WkJ<J>(/k 3 )<r(/ k2 )<l>*(/ k4 )> 

= (2n)- 2D+2d J d Tl dT 2 J dPdP'E[6 (y - x - b (n)) 5 (y' - x' - b' (r 2 )) 
exp (-£(P - kl ) 2 - £(P - k 2 ) 2 - £(P' - k 3 ) 2 - £(P' - k 4 ) 2 ) 
exp ( - I J ri J T1 P^PvPpD^P (b (s) - b (*')) dsds' 
-3 Jo' 2 Jo' 2 P'^PlP'^o (b' (*) - b' (s')) <W 

-3 Jo' 1 Jo' 2 P»PvP' p P'„D^°r (x - x' + b (5) - b' (*')) dsds')] + (1, 2 ^ 3, 4) 

(28) 

where the last term means the same expression with exchanged wave numbers. 
We introduce the spherical coordinates on the (n, r 2 )-plane t\ — rcosO and 
r 2 = rsin#. Let us rescale the momenta k = p^/r, k' = p'y/r , K = Pr2~i 
and K' = P'r^~i . Then, we can see that the four-point function (27) takes the 
form 

($(x)$(x')$*(y)<I>*(2/)) 

= / dedrrr- d -^-^ D -^F 4 (6, r~3 (x - y), r~i (x' - y'), 
r~5( x ' - x),r~2(x' - y),r~3(y' - x), 

r-i+i(X- Y),r-i+i(X'- Y'),r-i + i(X'- Y),r-i+i(X- Y')) 
It follows just by scaling of coordinates (the r-integral scales as twice the r- 
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integral in eq.(18)) that at x = x' = y = y' = the correlations are scale 
invariant with the same scaling dimension as in eq.(20). 

It is clear from eq.(28) that in the same way as we did it in eqs.(24)-(25) we 
can obtain finite upper and lower bounds on the correlation functions (28) by 
means of the Jensen inequalities. 

We could continue with higher order correlation functions. Again through 
an introduction of spherical coordinates in the (ti,...,ts) space we can show 
that 



scales with the same dimension as in eq.(20). The scaling of higher order corre- 
lation functions is now evident. We introduce the spherical coordinates for the 
r-integrals. The resulting scaling is a consequence of the fact that the r-volumc 
and P integrals have the scaling dimensions proportional to the order of the 
correlation function. 

4 A projection to d dimensions 

There is still another option that we let all X = Y = 0. In such a case 



By a scaling of momenta we can bring the propagator of eq.(3) to the form 



$(s3)**(yi)....**(tte)> 



(30) 



<# T (x,y)) = (2n)- D + d J dP 
E[6 (y-x- y/7b(l)) 

exp (-\t 2 - 2 ^ jl P»P a P v P p D^°<> (b (a) - b («')) dsds')] 



(31) 




(32) 
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Hence 

(A- 1 fry)) = J R|x-y|- d + 2 -^- d )( 1 -^ (33) 

where R is a positive constant. Hence, if all the correlation functions are scale 
invariant then 

d-2 . (D — d)(l — -v) 

$(x, 0) ~ A— + s $(Ax, 0) (34) 

We can prove all the inequalities of sec. 3 in this model. So, the upper bound 

for the two-point function reads 
IM-Hx,y))| = 

< 2 J °° dr J Q ds J S da' / du 1 du 2 dPT~ip(s' , ui)p(a - a', u 2 - ui) (35) 
P (l - a, r-i (y - x) - u 2 ) cxp (-i^lp^P a P v P p D^P ( Ul - u 2 )) 

The lower bound takes the form 
IM-Hx.y))! > (2n)- D + d f Q °°dTT- i JdP 

exp ( - ±t 2 -^ Jo P^P a P u P p E[D^P (a (r"5y - r"5x, a) - a (r"5y - r"5x, a') )dada']) 

(36) 

where the expectation value in the exponential on the r.h.s. of eq.(36) is equal 
to 

/ duj ds J* ds' (2™ (s,s'))-l exp (- .^byu 2 ) (3?) 
|u — t~2 S (y — x) + r _ 2 a' (y — x) |~ 47 

where w(a, a') = (a — a')(l — a + a'). It is finite if 7 < |. The bounds (35)-(36) 
in fact have the form 

i?i < |x-y| d - 2 +( I3 - d )( 1 - 7 )(^- 1 (x,y)) <i? 2 (38) 

with certain positive R\ and i? 2 . 
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The inequalities (38) can be proved from the inequlities (35)-(36) just by 
rescaling of variables. It is more tedious to show that the constants R\ and R2 
are finite and not zero (but the estimates reduce to finite dimensional integrals 
and are straightforward). 

We can project now to R d higher order correlation functions 

($(x)$(x')$*(y)$*(y')> 

= (^l- 1 (x,y)^- 1 (x',y')> + (x^x') 

= (2K)- 2D+2d J d Tl dT 2 J dPdP' 

£[5(y-x-b(r 1 ))5(y'-x'-b / (r 2 )) 

exp ( - ± J T1 J T1 P^P a P v P p D^P (b (a) - b (s')) dads' 

-I lo 2 Jo* P',P' a PlP' p D^ (b' (a) b' {*')) dsds' 

-3 Jo 1 lo 2 P»P»P' P P'«D^ (x - x' + b {a) b' (5')) <w)] + (a - x') 

(39) 

where (x — > x') means the same expression in which x is exchanged with x'. We 

introduce the spherical coordinates on the (n , T2)-plane n = r cos 0, t 2 = r sin 6* 

and we rescale the momenta k = p^/r, k' = p'y/r , K = Pri~2 . Then, we can 

see that the four-point function (39) takes the form 
($(x)$(x')$*(y)$*(y')> 

= )^rrr- d -( 1 -T)( D - d ) J F 1 4 (6',r-^(x-y),r-5(x' - y'), (40) 
r~5(x' - x),r~i(x' - y),r~i(y' - x )) 

The upper bound now reads 
|<$(x)$(x')$*(y)$*(y')>| < 

2(2ir)- 2D+2d J d Tl dT 2 Jo ds f° ds' / dPdP'£[<5 (y - x - V^b (1)) 
5 (y' - x' - ^b' (1)) 

exp ( - T ±P^P a P„P p D^P (Vnb (a) - Vnb (s')) (41) 
-llp'^P^D^P (V^b' («) - V^b' (a')) 

- T -^P p P v P' p P' a B^P (x - x' + Vnb (s) - V^ b ' )] + (x -» x') 
The expectation value (41) can be expressed by the transition functions (as usual 
for the Wiener process). The bound is scale invariant and the scale invariant 
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function on the r.h.s. could be calculated explicitly. The lower bound takes the 
form 

($(x)$(x')$*(y)$*(y')) 

> {2n)- 2D+2d J d Tl dT 2 J dPdP' 

cxp ( - E[\ J T1 J T1 P,P a P v P p D^P (a {a) - a («')) dsds' (42) 

\ J" Jo 2 PlKKP' P D ^ ap ( a ' M a ' (*')) d«fc' 

+ 5 Jo" Jo" P,P,P' p P' a D^P (a ( S ) - a' (,')) d^']) + (x - x') 

Here 

S S 

a(s) = x + (y - x)— + Vnc(- ) 

n ti 

and 

a'( S )=x' + (y'-x')- + ^c'(-) 

In the exponential of the formula (42) we have the expectation over three Gaus- 
sian processes. The first is with the mean (y — x)(s — s')/t\ and the covariance 
(s — s')(ti — s + s')/t\, the second has the mean (y' — x')(s — s')/t2 and the 
covariance (s — s')(r 2 — s + s')/t2 , the third has the mean x — x' + (y — x)s/ti — 
(y' — x')s'/t 2 and the covariance s(ti — s)/ti + s'(t 2 — s')/t 2 . The lower bound 
can be explicitly calculated and is given by a scale invariant function. It is 
clear how to calculate the higher order scale invariant functions and their scale 
invariant lower and upper bounds. 

5 Discussion 

The model discussed in sec. 3 is invariant under Euclidean rotations in D — d 
dimensions. Euclidean fields with the Osterwalder-Schrader positivity cannot 
be more regular than the free field (this follows from the Kallen-Lehman rep- 
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reservation) . In D dimensions the short distance behaviour of the correlation 
functions (18) is more regular than the one of the free fields. However, after 
setting all x = the behaviour is more singular than the canonical one in D — d 
dimensions. We can suggest a lattice model whose formal continuum limit co- 
incides with our scale invariant Euclidean field theory. The simplest possibility 
is to take the conformally flat metric placed on a sublattice just between the 
lattice sites of the scalar field (as the gauge fields in ref.|lO|). It seems that 
the Gaussian model of sec. 3 mixing the ferromagnetic and antiferromagnetic 
couplings would fail to be reflection positive (in any case this would not be easy 
to prove, see ref . ]Tl[| ) . The model (11) with the metric which is a square of a 
Gaussian field can be reflection positive (it is scale invariant with correlation 
functions expressed by the characteristic function J^(h) (13)). The continuum 
limit and the subsequent analytic continuation to Minkowski space would give 
a model of relativistic quantum field theory satisfying all Wightman axioms. 
The Wick square of a Gaussian field is an example of an infinitely divisible field 
| fl2| . An infinitely divisible field can take non-negative values. Its characteris- 
tic function has an explicit integral representation. Such a field can be a good 
candidate for a random metric. 

More interesting are the models in sec. 4. The lattice version of the La- 
grangian 

L = g^M- — — $* + V X $V X $* 

y v ' dX^ dX v 

will have the form —L = F + OF + MOM where is the reflection in the 
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plane perpendicular to one of the coordinates x (which will be chosen as time). 
This representation holds true because the random metric does not mix the 
temporal coordinates in V X $V X $*. Then, the reflection positivity results (see 
[ pr| ). In the lattice approximation we have to replace the (formal) Gaussian 
measure with a negative definite covariance (16) by a non-Gaussian measure on 
the metrics which has a formal Gaussian limit (e.g., replacing ^x 2 by 1 — cos a;). 
There remains to be proven that such a lattice approximation is convergent to 
the continuum. 
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